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Modeling Global Macroclimatic Constraints on 
Ectotherm Energy Budgets 1 

Bruce W. Grant 2 and Warren P. Porter 
Department of Zoology, University of Wisconsin, Madison, Wisconsin 53706 

Synopsis. We describe a mechanistic individual-based model of how global macroclimatic 
constraints affect the energy budgets of ectothermic animals. The model uses macroclimatic 
and biophysical characters of the habitat and organism and tenets of heat transfer theory 
to calculate hourly temperature availabilities over a year. Data on the temperature depen- 
dence of activity rate, metabolism, food consumption and food processing capacity are used 
to estimate the net rate of resource assimilation which is then integrated over time. We 
detail and explore the significance of assumptions used in these calculations. 

We present a new test of this model in which we show that the predicted energy budget 
sizes for 1 1 populations of the lizard Sceloporus undulatus are in close agreement with 
observed results from previous field studies. This demonstrates that model tests are feasible 
and the results are reasonable. Further, since the model represents an upper bound to the 
size of the energy budget, observed residual deviations form explicit predictions about the 
effects of environmental constraints on the bioenergetics of the study lizards within each 
site that may be tested by future field and laboratory studies. 

Three major new improvements to our modeling are discussed. We present a means to 
estimate microclimate thermal heterogeneity more realistically and include its effects on 
field rates of individual activity and food consumption. Second, we describe an improved 
model of digestive function involving batch processing of consumed food. Third, we show 
how optimality methods (specifically the methods of stochastic dynamic programming) may 
be included to model the fitness consequences of energy allocation decisions subject to food 
consumption and processing constraints which are predicted from the microclimate and 
physiological modeling. 

Individual-based models that incorporate macroclimatic constraints on individual resource 
acquisition, assimilation and allocation can provide insights into theoretical investigations 
about the evolution of life histories in variable environments as well as provide explicit 
predictions about individual, population and community level responses to global climate 
change. 



Introduction 

Impending global climate change will 
undoubtedly make the world less hospitable 
to some species of organisms and more hos- 
pitable to others. This realization has stim- 
ulated considerable interest in understand- 
ing how organisms respond to their changing 
environments and how these responses affect 
population size and dynamics. If we under- 
stand the mechanisms of interaction 
between environments and individuals, we 
can use individual-based models to predict 
the emergent properties (serisu Salt, 1979) 
of populations of interacting individuals and 
thereby model population level responses 
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to specific environmental changes (Huston 
et ai, 1988). Further, as changes in age 
structure and/or population size occur, we 
can also assess potentially destabilizing 
community-wide effects among other pop- 
ulations whose dynamics are ecologically 
interrelated. 

Predicting the outcomes of various sce- 
narios of climate change on populations of 
interacting individuals has traditionally been 
viewed as an applied area of ecology. How- 
ever, understanding the links between envi- 
ronments and individuals is of great theo- 
retical significance as well. Recently, several 
reviews have emphasized that investiga- 
tions of the mechanisms of individual 
response belong in the realm of more gen- 
eral questions about the evolution of life 
histories in variable environments (Cong- 
don et ai, 1982; Dunham and Miles, 1985; 
Dunham et ai, 1988a, 1989; Congdon, 
1989; Wilbur and Morin, 1988). A general 
approach to the question^of how to model 



154 

BEST AVAILABLE COPY 



Modeling Energy Budget Constraints 



155 



population response 10 climate change is to 
emphasize that selection acts on the 
responses of individuals to their environ- 
ments and ask how changing environmental 
constraints influence the expression of life 
history characters and the evolution of life 
history strategies. According to Dunham et 
al. (1989), the problem is to define the set 
of driving processes at the individual level 
which affect birth, death, emigration and 
immigration rates which are the system state 
variables and then predict the population 
size at the next time interval (see also Hus- 
ton et ai t 1988). The set of responses, in 
what Dunham et al (1989) refer to as a 
"time-ordered sequence of age-specific allo- 
cation decisions," by individuals deter- 
mines these state variables and is charac- 
teristic of the life history phenotype. Given 
any set of fixed or stochastically varying 
environmental constraints, and given the 
physiological state of the organism (which 
has resulted from past constraints and allo- 
cation decisions), the life history strategy 
transduces energy and mass into progeny. 
Heritable variation in phenotypic response, 
i.e.. life history strategy, that results in rel- 
atively greater ability to survive and repro- 
duce will be selected (see Congdon, 1989; 
Dunham et aL 1989). Therefore, the applied 
issue of modeling population responses to 
climate change is in fact a special case; of 
the general problem of modeling life history 
evolution in variable environments. 

Our aims in this paper are to describe how 
we model global macroclimatic constraints 
on animal eneifcy budgets, to address the 
implications and limitations of the assump- 
tions of our current modeling formulations, 
and to outline necessary data and tech- 
niques for more advanced models. We pre- 
sent a simple mechanistic model that estab- 
lishes bounds to the problem of predicting 
animal responses to variable environments 
for the case of small ectotherms (see Porter, 
1989, for a discussion of energy budgets for 
endotherms and larger ectotherms). We also 
examine concordance between model pre- 
dictions and observed geographical varia- 
tion in life history characters among 1 1 pop- 
ulations of Sceloporus undulatus to 
demonstrate that model tests are feasible 
and the results are reasonable. A compari- 



son such as this is relevant for predicting 
the effects of climate change on the life his- 
tories of small ectotherms since the mag- 
nitude of variation in environmental char- 
acters among sites differing in elevation by 
a few hundred meters or a few degrees in 
latitude is on the order of predicted varia- 
tion according to many scenarios of global 
climate change (COHMAP, 1988). 

An important goal of this paper is to make 
macroclimatic constraint/energy budget 
models more accessible to researchers in a 
wide variety of areas not presently involved 
with this type of modeling, and to provide 
sufficient direction for interested readers to 
construct their own versions of these mod- 
els. Additionally, we highlight the areas 
where our models are deficient and provide 
directions and guidelines for making spe- 
cific improvements tailored for use in other 
specialized implementations of the general 
questions we pose. 

Figure 1 provides a summary schematic 
of our modeling efforts and illustrates the 
major areas of emphasis in this paper to be 
detailed in the next section (Current Mod- 
eling Formulation). The first step is to model 
the mechanistic links between an animal and 
its environment and establish the con- 
straints that operate on its life history phe- 
notype. Any environmental factor that elic- 
its or constrains individual phenotypic 
responses and thereby affects birth and/or 
death processes is defined as operative and 
should be included (Spomer, 1973; Dun- 
ham et al t 1989). In this figure, we focus 
on modeling time constraints from the ther- 
mal biophysical environment on resource 
acquisition, assimilation and allocation; 
however, we suggest that including the effects 
of the resource, predation and social/demo- 
graphic environments is essential to mod- 
eling the complete system dynamics (see 
Dunham et al, 1989, for a more complete 
system diagram). 

The upper portion of Figure 1 indicates 
that macroclimatic input parameters (such 
as global solar radiation and local param- 
eters specific to the habitat of interest such 
as substrate type and topography) are cqm- 
bined with organism biophysical and phys- 
iological properties to determine the distri- 
bution of available thermal microclimates 
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(methods described below). This distribu- 
tion is defined as the operative biophysical 
environment, and all activities (such as 
thermoregulation, foraging, digestion, mate 
acquisition and predation) and inactivity are 
constrained to occur within this distribution 
(Grant and Dunham, 1988; Grant, 1990). 

The pattern of behavioral time allocation 
determines the daily accumulation of 
resources (energy, water and nutrients) as 
well as the daily resource expenditure due 
to maintenance and activity. Body temper- 
atures and activity times result from (a) 
tradeoffs among these resource gains and 
losses, (b) performance constraints stem- 



ming from temperature-dependent capaci- 
ties to digest food, locomote, and perform 
other vital functions (Huey, 1982), and (c) 
the constraint that thermally suitable 
microclimates are available in sufficient fre- 
quency in the ectothernTs home range for 
activity to occur (Christian et al t 1983; 
Waldschmidt and Tracy, 1983; Grant and 
Dunham, 1988, 1990; Grant, 1990). 

We define the total amount of assimilated 
resource available for allocation to growth 
(G), storage (S) or reproduction (R) over 
some "biologically meaningful time period" 
(Congdon et aL, 1982) in a life history as 
the discretionary resource. This amount is 
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similar to the "production budget" in Cong- 
don et at. (1982) except that we advocate 
the inclusion of reproductive activity 
expenditure in the fraction, R, that is allo- 
cated to reproduction. This inclusion is a 
necessary step in order to formulate indi- 
vidual reproductive effort according to 
< 1 9S °; see al *o Hirshfield and Tinkle, 
1975). In addition, for any temporal pattern 
of activity and T b selection, and for any allo- 
cation pattern of assimilated resource into 
G, S and R, there will be a combined asso- 
ciated nsk of mortality for any individual 
time and resource allocation phenotype 
(Dunham et al, 1989). This combined risk 
will be the product of all independent risks 
associated with each time or resource allo- 
cation "decision." 

The allocation of the discretionary 
resource among the competing categories of 
growth, storage and reproduction is a 
behavioral optimization problem and is 
hypothesized to maximize the expected sum 

?£fr I S" t ^ d future fecund ity (Williams, 
1 966; Hirshfield and Tinkle, 1 975). For sto- 
chastic or probabilistic environments, the 
problem of constrained optimization of 
resource allocation can be approached using 
standard optimization procedures (e.g., King 
and Roughgarden, 1982; Sibly et at, 1985* 

3 01 ?' \ 9 JL$ 5 om et al " ^90; Mangel and 
Clark, 1988; Taylor and Gabriel, 1992). 

Behavioral and physiological "choices" 
affecting the expected age-specific survival 
and fecundity of individual life history phe- 
notypes will determine population level 
shifts in age- and size-structure, generation 
time, population size and growth rate, and 
other emergent characteristics (Salt, 1979* 
Huston etal., 1988). These in turn will influ- 
ence the constraints that operate on the 
optimal life history phenotype stemming 
trom Other Environments" in Figure 1, 
e.g., by affecting the intensity of resource or 
mate competition or predation risk (Dun- 
ham et al % 1989). 



Current Modeling Formulation 
This section provides a brief overview of 
the methods we have developed to simulate 
operative environmental temperature, T 
availability and concomitant constraints ori 
the body temperatures, times of activity, 



and energy budgets of a ty pical cctotherm- 
a small rock dwelling lizard (see Porter, 
1989, for endotherms and larger ecto- 
therms). The most accurate method to assess 
T e availability and activity constraints is to 
measure these directly (e.g., Grant, 1990* 
Grant and Dunham, 1990). However, these 
methods are not useful in examining ther- 
mal constraints imposed by environments 
in the past or by hypothetical future envi- 
ronments under various scenarios of cli- 
mate change. Our objective is to bound the 
problem within the narrowest possible lim- 
its (Porter, 1989). This simulation proceeds 
in three parts: (A) delineate the physical 
microclimate characters for a hypothetical 
perching location, (B) calculate the opera- 
tive environmental temperature, T e , for a 
lizard with specified surface heat transfer 
properties, body posture and solar orien- 
tation, and (C) estimate the resource inges- 
tion and assimilation rate from regressions 
derived in lab as a function of body tem- 
perature, T b . 

Figure 2 illustrates the topographic com- 
plexity and concomitant thermal heteroge- 
neity of real environments (upper panel) 
which should be compared with the 
abstracted version of this environment 
required by this simulation (lower panel). 
This schematic illustrates what is needed to 
solve the energy balance equation for the 
substrate and for the hypothetical lizard. 
Simplifying assumptions in the microcli- 
mate simulation are adequate to reduce the 
complexity of the real environment into a 
small number of characteristic parameters 
and mechanisms of heat flux according to 
basic transient heat transfer theory (see also 
Kreith, 1973; Porter, 1989). Figure 3 illus- 
trates the energy balance equation which is 
used both for microclimate and ectotherm 
temperature simulations with only minor 
modifications. (see Porter et al., 1973, for an 
expanded and detailed treatment of these 
equations). For example, for the substrate, 
the diagonal heat balance equation simply 
excludes the term Tor metabolism, and for 
points below the substrate surface, all terms 
in the diagonal equation would be equal to 
zero except for conduction and storage of 
heat. The resulting simulated biophysical 
environment includes the range of available 
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surface temperatures for activity as well as 

£v« r^ rat °£ inacti ve lizards in deep 
caves. This will sufficiently bound the prob- 
lem by placing upper and lower limits on 
activity time, resource acquisition and 
assimilation. 



Physical microclimate simulation 

Real lizards can be active on horizontal 
or sloping rock surfaces and exposed to direct 
and/or indirect sunlight, partial shade, or 
deep shade exposed to the sky; alternatively, 
the lizard may cease activity and retreat into 
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Fig. 3. Overview of the coupled equations for the heat balance, along the diagonal, and for the mass balances 
along the horizontal t for a typical ectotherm (see Porter, 1989, for a more detailed treatment). Food ingested (at 
left) is processed and then available for allocation to growth (G), storage (S), or reproduction (R) in an ecological 
energy and water budget (at right): ■ 



caves and crevices to find T e s anywhere 
between the air temperature at 2 meters and 
the temperature in deep caves (Porter et ai t 
1973; Porter, 1989). However, this model 
only considers the bounding values of avail- 
able T e by simulating perch microclimate 
conditions in one of three extremes— a hor- 
izontal perch, a perch in full shade but 
exposed to the sky, or a perch in a deep cave 
(Fig. 2, lower panel). 

The simulation proceeds by specifying 
physical characters of the locality to be 
modeled such as elevation, site slope angle, 
substrate solar absorptivity, conductivity, 
and surface roughness (see Porter, 1989, for 
representative values). Maximum and min- 
imum air temperatures (TJ and wind speeds 
(V) at reference height (2 m), as well as the 
time of day when these extremes occur, are 
obtained from standard climate data sets for 
a particular date of interest. These inputs 
are fit to sinusoidal curves to obtain hourly 
estimates of T.(t, 2 m) and V(t, 2 m). The 
surface roughness parameter is used with 
the V(t, 2 m) curve to fit a logarithmic func- 
tion to estimate the wind speed V(t, z) as a 
function of height (z) and time (t) (Porter et 
al, 1973). To obtain T,(t, z) the simulation 
assumes the deep cave temperature is fixed 



at mean monthly temperature and estimates 
the surface temperature based on numeri- 
cally solving a version of the heat balance 
equation (Fig. 3) for the substrate and points 
below. This solution relies on a subroutine 
(SOLRAD, McCullough and Porter, 1971) 
to calculate the incident diffuse and direct 
global solar radiation load for the particular 
perch slope and azimuth angles. 

With the reference height and deep'cave 
temperatures known, the simulation uses a 
finite difference heat transfer algorithm 
(Adams Predictor-Corrector numerical 
integrator, G.E.A.R. package, Argonne 
Laboratories) to find the surface and inter- 
mediate temperatures between the reference 
height and "deep soil temperatures, T a (t, z) 
and T,(t, z). This simulation must be 
repeated for each of the two exposed perches 
(horizontal perch and a perch in shade) for 
which the T.(t, z) curve for the latter is cal- 
culated with the direct component of solar 
radiation from SOLRAt) set to zero. With 
the microclimatic conditions of radiation, 
temperature, and wind speed calculated for 
the hypothetical perches, we are ready to 
place an hypothetical lizard in these perches 
and estimate the operative environmental 
temperature, T e . 
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Fig. 4. The predicted range of environmental tem- 
peratures (T„max, T ei min, and T tf cave) available to 
active and inactive lizards for each hour of the day. 
Lizard activity is assumed to occur if either of the 
exposed microclimate T f s are within the user-defined 
range of T^max and T^min. During activity the lizard 
is assumed to regulate its T» (solid points) as close as 
possible to the "^preferred (see text for method of 
estimation). During inactivity, the lizard's T b equals 
the T„cave. 

Ectotherm operative environmental 
temperature simulation 

Surface heat-transfer properties and tem- 
perature dependent metabolic and water loss 
rates of the lizard of interest are used in 
combination with the driving microclimatic 
conditions to solve the diagonal heat bal- 
ance equation in Figure 3. This calculation 
is repeated at hourly intervals for each of 
the three perches in (A) which will delineate 
the range of available T e s to active and inac- 
tive lizards (Fig. 4). Since the exposed hor- 
izontal perch represents the maximal T c 
available, the lizard is assumed to maintain 
the greatest possible fraction of its body sur- 
face exposed to sunlight throughout the day. 
Body orientation is assumed to have no 
effect on the T e s for the exposed shaded and 
deep cave perches since there is no direct 
sunlight. The T e for the deep cave perch, 
which will be used during inactivity, is 
assumed to equal the deep cave tempera- 
ture. 

Upper and lower acceptable bounds for 
body temperature, T b , during activity are 
input for the lizard of interest, and for each 
hour of the day the lizard will be inactive 
if both exposed T e (t)s are outside the range 
of acceptable T b (both too hot or both too 
cold), otherwise the lizard is active (see Fig. 
4). With bounds set for both the length of 



time and for the range of body temperatures 
during activity, we are ready to allow the 
hypothetical lizard to attempt to select its 
physiologically optimal body temperature 
and consume and digest food resources. 

Ectotherm physiological simulation 

Laboratory determined temperature- 
dependent metabolic rate, food consump- 
tion and digestive processing rate are used 
in combination with activity and T b esti- 
mates from (B) to estimate the total amount 
of energy (less maintenance and digestive 
costs) ingested and assimilated by the lizard. 
Resting metabolic rates for this simulation 
are obtained from a regression in Bennett 
and Dawson (1976; see Appendix A). 
Although this equation was calculated from 
a wide variety of lizards and represented an 
important early contribution to studies of 
physiological allometry, the predicted value 
as a function of T b and body size might differ 
considerably from the actual measured value 
for the particular lizard of interest due to a 
wide variety of important ecological and 
evolutionary considerations (Bennett and 
Dawson, 1976; McNab, 1979; Nagy, 1982). 
Ideally, the equation relating metabolic rate 
to T b and body size should be estimated 
from data on individuals from the Specific 
population and possibly also season of 
interest. The metabolic rate of an active liz- 
ard is assumed to be a constant fraction 
times the resting rate (s activity scope, Ben- 
nett and Dawson, 1976; Congdon et ai t 
1982; see Appendix A). 

Maximal food consumption and assimi- 
lation rates were determined for a lizard as 
a function of T b by estimating the steady 
state rates of ingestion and excretion (J/day) 
by lizards held at constant cage temperature 
in an ad lib feeding regime (data from See- 
toporus undulatus, Waldschmidt, 1984; 
Waldschmidt et al, 1987). Consumption 
and assimilation rates were fit to polyno- 
mial regressions, and the energetic cost of 
digestion (referred to as the specific dynamic 
action) was estimated in lab as a function 
of T b (see Appendices A and B). To simulate 
the consumption and assimilation rates for 
the hypothetical lizard, it is assumed that 
these rates for an hour of activity at a given 
T b equal V 24 of the laboratory estimated daily 
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rates at that T b . This is equivalent to assum- 
ing that these lizards are entirely "process 
limited 11 (Congdon, 1989), i.e., their energy 
budgets are not constrained by their inges- 
tion rate, but rather by their digestion and 
assimilation rate. Such an assumption would 
represent an upper bound on the amount of 
energy consumed. However, since this 
assumption may not hold for lowland desert 
environments^., Dunham, 1981; Cong- 
don, 1989), the simulation allows for inges- 
tion limitations by simply taking a constant 
fraction of the maximal ingestion rate 
obtained from the ad lib feeding regime such 
as 75% or 50% of maximal food consump- 
tion capacity. 

Perhaps a more significant assumption of 
the present physiological model is that no 
feeding is allowed to occur when the animal 
is inactive, and thus no digestion and assim- 
ilation can occur either. The fraction of food 
eaten during any particular hour of activity 
is assumed to be assimilated during that 
hour according to the regression for assim- 
ilation rate as a f(Tb). In other words, regard- 
less of what T b the lizard exhibits in sub- 
sequent time intervals during which 
digestion would normally occur in real liz- 
ards, the mass budget of the hypothetical 
lizard is credited with digesting the con- 
sumed food, and debited the energetic cost 
of digesting that food, depending on its T b 
during the time interval of consumption, 
i.e. t assimilation is instantaneous. This could 
represent an upper bound to the processing 
capacity of the lizard, since in miany situa- 
tions the temperature during inactivity 
(overnight or at midday) is lower and assim- 
ilation rate would be greatly reduced (see 
Appendix A). On the other hand, this model 
could underestimate food processing capac- 
ity for lizards that feed in cool microcli- 
mates such as in early morning (e.g., desert 
sceloporines) or underwater (e.g., marine 
iguanas) and have access to warmer micro- 
climates later in the day in which to digest 
their food at a thermally enhanced process- 
ing rate. This assumption is a current focus 
of concern and in a later section we describe 
how to allow for batch processing of ingested 
resource in which the lizard may digest and 
assimilate food during time intervals sub- 
sequent to consumption as a function of its 
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Fig. 5. Sample output from the physiological model 
for a 5 g Sceloporus undulatus describing the temper- 
ature dependence of resource intake, processing, and 
expenditure (in J/hr, equations in Appendix A). 
Qavaiuioo« and Qavaiujo* represent the total assimilated 
energy for 100% and 50% of ad lib. feeding, respec- 
tively, Qmet is the total metabolic expenditure for an 
activity scope of 1.5 times resting. The difference 
between Q MCT and either Q AVA il is the net amount of 
assimilated energy available for allocation in the energy 
budget. Arrows indicate the T b at which net assimilated 
energy is maximal for each feeding regime, and the 
model uses this temperature as the preferred temper- 
ature during activity by the hypothetical lizard (sec Fig. 



T b at these later times. Batch processing 
models represent an important step in 
incorporating thermal constraints on tran- 
sient resource processing capacity (see 
below). 

Figure 5 is a sample output from the phys- 
iological model for a 5 g Sceloporus un- 
dulatus and describes the temperature 
dependence of resource intake and expen- 
diture (J/hr, equations in Appendix A). The 
convex decelerating curves represent the 
total assimilated energy (consumption frac- 
tion assimilated) for unlimited food and for 
a reduced food level to 50% of ad lib., 
(Qavaiu too% and Qavaiu5o%» respectively). The 
exponentially increasing curve is the met- 
abolic expenditure, Q MET , for an activity 
scope of 1.5 times resting. The difference 
between the Q MET and either of the total 
assimilation curves represents the net 
amount of assimilated energy available for 
discretionary allocation to growth, storage, 
or reproduction for either of the feeding 
regimes. As can be seen, these data indicate 
that an hypothetical reduction in available 
food from ad lib. conditions to only 50% of 
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ad lib. food consumption greatly reduces 
discretionary energy at any body tempera- 
ture. Note that at temperatures above about 
36°C for the 50% ad lib. regime, the lizard 
would be in a negative energy balance. Dur- 
ing activity, the hypothetical lizard is 
assumed to have access to an intermediate 
thermal microclimate in which it may be 
active at the T b which maximizes its net 
energy discretionary energy intake per hour 
(see "T b ,preferred" in Fig. 4, and the arrows 
in Fig. 5). Note in Figure 5 that this tem- 
perature is lower when the lizard has limited 
access to food than when food is unlimited. 
This model is consistent with results from 
other ectotherms such as freshwater fishes 
(Brett, 1971;Magnuson and Beitinger, 1978; 
Crowder and Magnuson, 1983) and aquatic 
insects (Vannote and Sweeney, 1980), but 
has yet to be tested in lizards. 

Model Predictions 
Development and testing of the micro- 
climate and steady-state ectotherm temper- 
ature simulations have been ongoing for the 
past two decades. These studies have shown 
how coupled daily and seasonal biophysical 
constraints on activity time and resource 
acquisition and processing capacities affect 
(a) individual daily thermoregulatory 
behavior (Beckman et ai. t 1973; Porter et 
al. t 1973), (b) seasonal potential for growth, 
storage and reproduction (Porter and James, 
1979; James and Porter, 1979; Porter and 
Tracy, 1983; Porter, 1989), (c) species dis- 
tribution limits (Porter and Tracy, 1983), 
and (d) trophic interactions among species 
(Porters aL, 1975). In this section, we pre- 
sent an additional test of this modeling in 
which we compare observed and predicted 
integrated annual energy budget sizes for 
Sceloporus undulatus. which is one of the 
most well-studied lizards in North America. 

The extensive geographic range of S. 
undulatus (from Florida to New Jersey, west 
to Utah, and south into north-central Mex- 
ico) as well as their highly variable life his- 
tory characters among populations (see 
reviews in Tinkle and Dunham, 1 986; Dun- 
ham et ai t 1988a) render this lizard ideal 
for modeling macroclimatic constraints 
along environmental gradients. According 
to extensive analyses by Tinkle and Dun- 



ham (1986) and Dunham et aL ( 1 988a) using 
data from 1 1 populations throughout this 
range, very little of the observed life history 
variation can be explained by any simple 
hypothesis (due either to absence of pattern 
or absence of data). In particular, these 
authors could find no evidence to support 
the suggestion by Ferguson et aL ( 1 980) and 
Ballinger et aL (198 1) that life history char- 
acters should show convergence in similar 
environmental types, which suggests that 
environmental constraints on these lizards 
exert little if any effect. However, this could 
easily have been due to the coarseness of 
the initial environmental classification sys- 
tem (Tinkle and Dunham, 1986), which 
relied heavily, on structural differences 
among habitats and would necessarily only 
provide a crude estimate of the daily and 
seasonal thermal environmental constraints 
impinging on these lizards. Direct estimates 
of variation in temperature availability are 
necessaiy to delineate macroclimatic con- 
straints on daily and seasonal activity time 
and concomitant resource accumulation and 
allocation in a life history (Grant and Dun- 
ham, 1988, 1990). 

The macroclimatic constraint/energy 
budget model we presented in the previous 
section provides a greatly improved test of 
the hypothesis that environmental factors 
affect variation in life historical energy bud- 
gets among populations of S. undulatus. 
Table 1 summarizes available life history 
data that we used to estimate differences in 
annual reproductive output among 10 pop- 
ulations compiled by Dunham et aL ( 1 988a) 
and one population in New- Jersey (P. Nie- 
wiarowski, personal communication). If 
annual net storage were near zero (which 
would be the case if lipids cycled but showed 
no annual accumulation), and if growth were 
near zero (which is reasonable for repro- 
ductively active adult lizards, Andrews, 
1 982), then this estimate of the annual mass 
allocation to reproduction for each site may 
be quite similar to the total annual energy 
budget size which is the final output pre- 
diction of our model. 

To drive the microclimate and ectotherm 
models for each of the 1 1 study sites, we 
used long term average air temperatures 
from nearby weather stations, we estimated 
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local substrate physical characteristics from 
the descriptions in the respective sources, 
and we used biophysical and physiological 
data on Sceloporus undulatus as described 
above (see Appendices A and B). Addition- 
ally, since body size greatly affects stomach 
capacity, food consumption and processing 
rate (Appendix A), we used the observed 
average snout- vent lengths of adult females 
for each site for these simulations. Lastly, 
it should be noted that we assumed, due to 
an absence of evidence, that these biophys- 
ical and physiological characters for S. 
undulatus do not differ among seasons and 
sites (we return to this assumption below). 

Figure 6 shows the observed annual total 
dry wt of eggs vs. the predicted maximum 
annual mass available for allocation to eggs 
according to our simulation. Note that none 
of the observed values greatly exceed the 
maximum predicted values for each site 
which shows that the model does not under- 
estimate the upper bound for any of these 
sites. Secondly, note that there is a signifi- 
cant polynomial regression of observed 
annual reproductive output, Y, on predicted 
annual energy budget size, X (Y = -1.20 
+ 1.87 X - 0.17-X 2 ,i?* = 0.842, n- 11, 
P < 0.01). Given the large number of 
assumptions required to generate the pre- 
dicted values, as well as the diverse and 
often crude field methods used to estimate 
the observed reproductive parameters 
(especially clutch frequency), some may find 
it remarkable that any relation exists at all, 
let alone one that explains this much of the 
variance. 

The comparison in Figure 6 suggests that 
environmental constraints among these sites 
do operate on daily and seasonal activity 
time and thereby exert an effect on observed 
annual reproductive output. This is broadly 
consistent with the earlier suggestion (see 
above) that similar life histories, at least with 
respect to total production budget sizes, may 
result in similar environments. However, 
this model can make no predictions about 
how this resource should be allocated in 
terms of egg size, clutch size, and the num- 
ber of clutches per year— all of which vary 
greatly among sites (Table 1) and are basic 
characters of a life history (Dunham et ai t 
1989). Our result only suggests that the total 
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Predicted Energy Budget Mass (gm/yr) 

Fic. 6. Observed annual reproductive output (g dry 
wt) vs. predicted annual energy budget size (g dry wt) 
for 11 populations of Sceloporus undulatus from 
throughout their geographic range. The "observed" data 
and habitat classes for all sites except New Jersey were 
summarized in Dunham et al. (1988a. see Table I). 
Data from New Jersey are from ongoing work by P. 
Niewiarowski (personal communication). Squares are 
from canyonlands (Left to right— Pinos Altos [New 
Mexico], Utah, Colorado, and Arizona), circles are from 
grasslands (left to right— Nebraska, Kansas, Texas and 
Lordsburg [New Mexico]), and triangles are from 
woodland sites (top to bottom— Ohio, South Carolina 
and New Jersey). A reference line appears with slope 
1,0 and intercept 0. The polynomial regression 
(Observed - -1.20 + 1.87 Predicted - 0.17- 
Predicted 2 , R 1 - 0.842) is based upon data from all 1 1 
sites. 

amount of resource available for allocation 
is constrained, which may provide insights 
into allocation rules when combined with 
other modes of inquiry (see Dunham et ai \ 
19886). 

An important function of broad compar- 
isons such as in Figure 6 is to direct future 
efforts by suggesting specific hypotheses to 
be tested by experiments (Dunham et ai t 
1 9886, 1989), For example, several hypoth- 
eses may be posed to explain why some of 
the observed values are much less than the 
maximum predicted and why the regression 
is a non-linear function. First, the model 
assumes that all of these lizards are limited 
by their capacity to process food with no 
limitations on ingestion rate (see Congdon, 
1989); consequently, greater "annual accu- 
mulation is predicted in warmer environ- 
ments with higher processing rate and lon- 
ger integrated annual activity time. 
However, if the observed warmer environ- 
ments also tend to have lower availabilities 
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Table I . Summary of available life history data on Sceloporus undulaius used to estimate observed geographic 
variation in annual reproductive output. 



Study site 


Habitat 
type 


Dutch size 


Dutches 
per year 


Eu wet mass 


Adult 
9SVL 


Weather data from 


Arizona* 


C 


8.3 


3 


0.29 


65 


Phoenix, AZ 


Utah" 


c 


6/3 


3 


0.36 


69 


Cedar City, UT 


Colorado* 


c 


7.9 


2 


0.42 


70 


Grand Junction, CO 


Pinos Alios, NM* 


c 


7.2 


2-3 


0.29 


63 


Beaverhead, NM 


Lordsburg, NM d 


G 


9.9 


4 


0.24 


68 


Lordsburg, NM 


Kansas' 


G 


7.0 


1-2 


0.26 


57 


Topeka. KS 


Nebraska' 


G 


5.5 


2 


0.23 


55 


North Platte, NE 


Texas* 


G 


9.5 


3 


0.22 


57 


San Angelo, TX 


Ohio* 


W 


11.8 


2 


0.35 


75 


Lancaster, OH 


South Carolina* 


W 


7.4 


3 


0.33 


63 


Aiken, SC 


New Jersey* 


W 


9.0 


2 


0.395 


73 


Pemberton, NJ 



* Modified from Dunham et ai t 1988a. Sources: • Tinkle and Dunham (1986). * Tinkle (1972). « Tinkle and 
Ballinger (1972).- Vinegar ( 1 9 7 5a. b). * Ferguson et ai ( 1 980). ' Ballinger efa/.(l981).»P. Nie wiarowski (personal 
communication). Habitat types are C-canyonlands, G-grasslands, W- woodlands. Weather data were obtained 
from the U.S. Climatological database for the recording station nearest and most similar in elevation to the 
exact location of each study site. Annual reproductive output was calculated as clutch sizeclutches/yregg wet 
mass, for which we used values of 1.5, and 2.5 for dutches/yr of "1-2" and "2-3", respectively. We used data 
in Vitt (1978) to convert to dry weight for comparison with the model prediction (dry weight [g] « 0.4743 -wet 
weight [g)). 



of food (such as may occur in desert grass- 
lands or canyonlands), then the observed 
values for energy budgets in these habitats 
will be much less than predicted, whereas 
in cooler food-rich environments (such as 
northern and montane grasslands, or wood- 
lands) observed and predicted values are 
expected to be much closer. These effects 
could lead to increased deviation of the 
observed reproductive output from the pre- 
dicted maximal energy budget size as the 
latter increases. In fact, the regression after 
omitting the wannest and driest site, Ari- 
zona, is highly linear (Y = - 0.052 + 0.912- 
X, R} = 0.868, P < 0.0001) with a slope 
not different from 1.0 (/ - 0.70, df - 8, P 
> 0.05). This suggests that the energy bud- 
gets of S. undulatus are limited by food 
ingestion rate at the Arizona locality, 
whereas, at the other sites, budgets are lim- 
ited by food processing capacity (see also 
Congdon, 1989). These results clearly point 
to the need to estimate food availability in 
field studies of life history variation (e.g., 
Dunham, 1978), and to conduct resource 
manipulation experiments to determine in 
which environments do increases in food 
availability actually lead to increases in the 
amount of resource available for allocation 
(e.g., Waldschmidt, 1983; Guyer, 1988a, b\ 
Grant and Contreras-Arquieta, unpub- 
lished data). 



There are several other hypotheses that 
can account for the non-linear regression in 
Figure 6, all of which underscore the need 
to incorporate other types of inquiries in 
developing and testing energy budget mod- 
els. For example, since warmer environ- 
ments may also be drier, such as in deserts, 
and since about 50% of the mass of freshly 
oviposited Sceloporus undulatus eggs is due; 
to water (Vitt, 1978), water may limit annual 
reproductive output more so than does 
energy in warm and dry environments. In 
other words, even if a desert lizard has suf- 
ficient energy to produce an additional egg, 
unless she also has the minimal water store 
to hydrate that egg she will have to allocate 
that energy elsewhere in her life history (e.g., 
to growth or storage). Modeling this sce- 
nario requires a coupled simulation of the 
annual water and energy budgets of the study 
lizards (we thank A. E. Dunham for this 
intriguing suggestion). 

Another explanation for the non-linear 
regression in Figure 6 is if there were a lower 
fractional allocation to reproduction at larger 
total energy budget size. This could result 
from a constraint on maximal egg produc- 
tion rate at larger potential production bud- 
get sizes (e.g., if the rate of yolking a clutch 
were limiting, or if egg development or post- 
hatching growth rates were low and ovi- 
positions were constrained to occur early in 
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the breeding season to insure juvenile sur- 
vivorship). On the other hand, if demo- 
graphic characters were such that low adult 
body sizes and high repoductive output were 
selected (e.g., high adult mortality and early 
age of maturity), concomitant total energy 
budget sizes would be low and the fraction 
allocated to reproduction could be high (dis- 
cussed in Tinkle and Dunham, 1986; Dun- 
ham et aL, 1988a). 

Some of the residual variation in Figure 
6 no doubt results from several basic sources 
some of which are peculiar to the model. 
For example, a problem with this simula- 
tion is that the driving environmental mac- 
roclimatic data were not available from any 
of the study sites either because these data 
were not collected or not reported. Typi- 
cally, the data that were available were from 
nearby airports (see Table 1), the biophys- 
ical environments of which may poorly 
reflect those available to lizards on their 
rocky and/or coarsely vegetated home 
ranges. Another potential source of varia- 
tion is that the model assumes clear skies 
throughout the year and undoubtedly the 
study sites differ in annual cloud cover. 
Cloud cover generally reduces environmen- 
tal temperatures which could decrease total 
energy budget size if these cooler tempera- 
tures inhibit activity, foraging and diges- 
tion; on the other hand, cloud cover could 
have the opposite effect on activity, food 
consumption, and energy budget sizes if 
sunny conditions are thermally stressful, 
such as during southern and low elevation 
summer days. We strongly suggest, as have 
others (Waldschmidt and Tracy, 1983; Tracy 
and Christian, 1986; Dunham et aL, 1989; 
Grant and Dunham, 1990), that studies of 
ectotherm population biology include 
detailed estimates of the thermal environ- 
ments available to and used by the study 
individuals. 

Another potential source of error in the 
model stems from the assumption of iden- 
tical thermal physiologies for the 1 1 5. 
undulatus populations modeled. What little 
data are available on activity T b suggest the 
opposite may be true. For example, S. undu- 
latus exhibit different field active average 
T b s in New Jersey than they do in Nebraska 
by several °C, and the effect of this on their 



activity and seasonal energy budgets may 
be great (P. Niewiarowski and W. Roosen- 
berg, unpublished data). However, com- 
parable field data are unavailable from other 
sites, and even simple comparisons under 
laboratory conditions are lacking to support 
or reject the use of the same metabolic and 
food processing capacity equations among 
all sites and seasons (Appendix A). Seasonal 
acclimation of physiology is well known for 
lizards (e.g., Nagy, 1 982), and when the gen- 
erating physiological regressions are known 
for each population and season, simulations 
of this type will be greatly improved. 

But a deeper source of variation may be 
rooted in unexplored aspects of the ecology 
and natural history of each system that 
directly affect specific variables in the model. 
For example, in populations with high den- 
sity or with a short and intense breeding 
season, reproductive activity expenditure 
can detract from the production budget 
(Anderson and Karasov, 1988) and the 1 1 
S. undulatus study populations differ greatly 
in density (see sources in Table 1). This 
argues for estimating season-specific activ- 
ity scopes (the fractional increase in meta- 
bolic rate above resting due to activity) for 
each site in order to correct the total main- 
tenance expenditure. As another example, 
if the predation environments vary among 
sites, constraints may be placed on activity 
and thermoregulation that could go either 
way depending on the nature of the preda- 
tor-prey interaction (e.g., high T b selection 
to maximize escape probability, or low T b 
due to selection of cool refuge microhabitats 
to reduce predator exposure). Additionally, 
high predation rates may select directly on 
energy allocation phehotypes by affecting age 
at maturity and age-specific reproductive 
effort (Stearns, 1976; Dunham etaL 1988a). 
Finally, there could be operative constraints 
due to phytogeny which could elicit com- 
plex associations among model parameters 
among sites, despite the lack of an overall 
effect of phylogeny in a simple test by Dun- 
ham et aL ( 1 988a). There is an obvious need 
for higher resolution in such analyses. 

Suggestions for Advanced Models 
In previous sections, we highlighted sev- 
eral areas where simplifying assumptions 
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used by the microclimate and physiological 
models deviated from biologically realistic 
levels of complexity. This section explores 
three such areas and presents ways in which 
the model presented above could be 
improved. We address how to incorporate 
biologically realistic levels of (A) thermal 
microclimate heterogeneity, (B) digestive 
performance, and (C) optimal allocation of 
time and resource to maximize fitness. 

Spatio-temporal thermal heterogeneity 
and activity constraints 

The microclimate model presented above 
uses an extremely simplified abstraction of 
the thermal environment to determine daily 
and seasonal constraints on activity and T b 
selection. Increased metabolic and feeding 
rates associated with activity are viewed as 
step functions varying between 0, for inac- 
tivity, and some maximum during activity, 
for which the edges of these steps corre- 
spond to variation in the ranges of available 
T e s (see Fig. 4). However, several previous 
researchers (e.g., Waldschmidt and Tracy, 
1983; Tracy and Christian, 1986; Grant, 
1990; Grant and Dunham, 198jB, 1990) have 
emphasized the importance of using the 
actual distribution of available thermal 
microclimates to delineate constraints on 
individual activity and home range use. 
Further, according to results in Grant ( 1 990), 
the actual shape of the T e (t) distribution may 
directly influence the cost/benefit econom- 
ics of lizard T b selection during activity. 

According to data from continuous obser- 
vations of free-ranging Sceloporus merriami 
from the Grapevine Hills iyi Big Bend 
National Park, Texas (Grant and Dunham, 
1988), individual activity rates (move- 
ments, feeding and social displays) during 
the morning activity period increased and 
decreased much more smoothly than do step 
functions. These authors concluded that the 
timing of peak activity rates as well as ran- 
dom microclimate use coincided with the 
time at which an indirect measure of the 
average available T e equalled the mean 
selected T b (calculated over all times in the 
morning). Subsequent work by Grant and 
Dunham (1990) greatly improved the res- 
olution of the thermal environment in this 
system (by using arrays of hollow-body 



painted models of S. merriami) and pro- 
vides a direct estimate of the distribution 
of available T e s for reanalysis here. 

In Figure 7 (upper panel) is plotted the 
fraction of model temperatures, T e , within 
the interval 26.0°-36.5°C which encloses 
95% of all T b data from active lizards on 
cloudless summer days in the Grapevine 
Hills, and is a measure of the range of 
acceptable temperatures for active lizards 
under these conditions at this site (modified 
from Grant and Dunham, 1 990). This figure 
indicates that by about 1 000 hr almost 100% 
of the habitat is thermally amicable to these 
lizards, whereas, before 900 hr and after 
1200 hr less than about 25% of the habitat 
is suitable for activity. Overlain on this fig- 
ure are the average rates of individual 
movements (solid symbols, modified from 
Grant and Dunham, 1988), from which it 
is apparent that the temporal pattern of nat- 
ural movement rates are poorly described 
by any step function as required by the 
model. Instead, movement activity shows 
more of a sinusoidal relation with time 
enclosing and in phase with the T e avail- 
ability curve. This suggests that the activity 
scope (sthe ratio of active to inactive met- 
abolic rate) should be a function of the dis- 
tribution of available T e within specified 
temperature bounds, i.e., when few micro- 
climates exhibit acceptable T e s for activity, 
lizard activity scopes will be low, whereas 
when most microclimate T e s are suitable for 
activity, scopes will be maximal. 

A similar result holds between the 
observed feeding rates and the distribution 
of T e within the interval 26.0°-36.5*C (Fig. 
7, lower panel), and somewhat surprisingly, 
the relationship is linear OR 2 - 6.841, P < 
0.01). This shows that individual feeding 
strike rates by these wait-ambush arthropod 
predators can be accurately predicted by 
estimating spatio-temporal thermal hetero- 
geneity. The explanation for the high degree 
of linearity in this relation may be as much, 
if not more, related to thermal environ- 
mental constraints on arthropod activity 
rates than to lizard thermoregulatory and 
prey capture strategies; however, data on 
this are lacking. 

The above results clearly suggest that any 
useful macroclimatic model to predict 
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activity and food consumption constraints, 
must at least attempt to simulate these con- 
straints at biologically realistic levels of 
complexity. One way to accomplish this, is 
to repeat the microclimate model numerous 
times for the same global site conditions but 
using randomly determined local perch and 
lizard body orientations. Each iteration 
would result in the T c (t) curve for one perch 
over a 24 hr period, and repeated iterations 
enable calculation of the distribution of 
available temperatures for the initial 
assumed topography (i.e., source distribu- 
tion for perch orientations). Figure 8 (upper 
panel) illustrates the results from 40 simu- 
lations using randomly determined perch 
locations and lizard orientations on a hemi- 
spherical rock in mid-summer. These data 
may be used to calculate the fraction of 
available microclimates and infer activity 
constraints as a function of time of day for 
any hypothetical activity temperature range 
(Fig. 8, lower panel). These calculations can 
be used to model time dependent activity 
scopes and foraging rates more realistically, 
and examine the impact of differences in 
patterns of movement rate and foraging 
mode on energy budget sizes among pop- 
ulations or species. 

Batch processing of ingested food 

Recent work by Penry and Jumars (1986, 
1 987) has demonstrated the utility of chem- 
ical reactor theory in modeling the reaction 
kinetics in animal guts. In this view, a gut 
can be modeled as a continuous flow through 
reactor (for animals that feed more or less 
continuously as with many filter feeders and 
herbivorous mammals), or as a batch reac- 
tor (for animals that feed in discrete pulses 
such as many wait-ambush carnivores). 
These models specify the; reaction kinetics, 
i.e., the rate equations relating ingested food 
to assimilable products, and simulate the 
steady state or transient mass balance for 
the gut. Variation in input parameters (e.g., 
morphology, physiology or feeding behav- 
ior) may provide inisights into gut design 
constraints and identify critical variables 
affecting optimal gut function (Sibly, 1981; 
Taghon, 1981; Taghon and Jumars, 1984; 
Trover, 1984; Penry and Jumars, 1987). 

Since ectothermic lizard activity and 
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Fio. 7. UPPER PANEL: The fraction of microcli- 
mate temperatures within the interval (26.0*-36.5*) 
available to Scetoporus merriami for a typical mid- 
summer day in the Grapevine Hills, Big Bend, Texas. 
This interval encloses 95% of all body temperatures 
for these lizards at this site. Superimposed on this plot 
are the reanalyzed data on individual movement rates 
from continuous observations of free-ranging S. mer- 
riami (solid points, Grant and Dunham, 1988). These 
data are depicted as **% of maximal rate" in order to 
use the same scale as the microclimate fraction data* 
LOWER PANEL: The number of feeding strikes per 
30 min by free-ranging S. merriami (from Grant and 
Dunham, 1988) is linearly correlated with the fraction 
of available microclimate . temperatures within the 
interval (26.0°-36.5 # ). Both panels indicate that tem- 
porally variable feeding rates and activity scopes need 
to be included in macroclimatic constraint/energy bud- 
get models. 

feeding rates are fairly discrete, depending 
on the temporal pattern of thermal micro- 
climate availability (e.g., see the previous 
section and Fig. 7), we suggest that a batch 
reactor model may provide an adequate 
representation of lizard gut functioning. Such 
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Fig. 8. UPPER PANEL: Simulated microclimate 
temperatures for 40 randomly oriented substrate sur- 
faces from repeated iterations of the microclimate pro- 
gram. Each orientation was randomly selected from a 
point on the surface of a fictitious hemispherical rock; 
however, any other source distribution, i.e.. rock shape 
may be used. LOWER PANEL: The fraction of sim- 
ulated microclimate temperatures a vailable to a hypo- 
thetical lizard within two typical body temperature 
intervals (26*-38% solid points, and 30»-36*. open 
points). We based these fractions on 120 calls to the 
microclimate program (the first 40 of which are shown 
in the upper panel). 

a model is needed to more realistically 
incorporate thermal constraints on resource 
processing capacity in any macroclimate 
constraint/energy budget simulation by 
allowing for digestion and assimilation to 
occur during times subsequent to consump- 
tion and at rates determined by the T b at 
these later times (whether or not the lizard 
is active). Thus, a hypothetical lizard that 
has recently emerged with a low T b may 
overfeed relative to its steady state ingestion 
rate at that T b , and process these resources 
later in the day. On the other hand, if inac- 
tivity T b ts too low for complete gut clear- 
ance during the inactive period, the lizard's 



consumption rate upon emergence for 
activity will be appropriately reduced. 

We have developed a simple model for 
batch processing of ingested food based on 
the general formulation in Penry and Jumars 
(1987) and using the specific digestive per- 
formance relations derived for Sceloporus 
undulatusztid presented in Appendix A. The 
model assumes that during each hour of 
activity the lizard (a) thermoregulates per- 
fectly (constant T b ), (b) ingests enough food 
to completely fill its gut (its energy budget 
is limited by its food processing capacity), 
(c) digests this food at an hourly rate of V 2 [ 
of the daily rate estimated in lab at the same 
T b (digestion rate is linear in time), and (d) 
metabolizes energy at a rate equal to the 
sum of its resting and activity rates as well 
as the cost of digestion. Also, during inac- 
tivity, the lizard does not feed and its hourly 
metabolic rate equals its resting rate plus 
the cost of digesting the food remaining in 
its gut, both of which are calculated for the 
inactivity cave T b . A detailed explanation 
of this simulation is presented in Appen- 
dix B. 

Results for several sensitivity analyses 
appear in Figure 9, each panel of which 
shows the cumulative energy budget size for 
the same set of default input parameters 
(upper curve) and the rate of energy accu- 
mulation due to a change in one of these 
parameters (lower curve). Note the decrease 
in the predicted energy budget size due to 
either lowering inactivity T b (Panel A) or 
activity T b (Panel B), These result from an 
increase in the food passage time and a con- 
comitant drop in the hourly rate of digestion 
associated with lower T b , and suggest that 
fed lizards should select warmer T b s. This 
is consistent with numerous lab studies 
demonstrating higher T b regulation after 
feeding by reptiles (reviewed in Wald- 
schmidt et al t 1987) as well as with field 
studies that suggest sleeping site thermal 
microclimate selection by lizards may be 
related to nighttime digestion requirements 
(Christian et ai t 1984; Huey et al. t 1989). 

Panel (C) illustrates the sensitivity of the 
energy budget size to an increase in activity 
rate (a 10% increase in activity scope pro- 
duced a 25% drop in the rate of energy accu- 
mulation). This difference in scope is prob- 
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Fio. 9. Predicted cumulative energy budget size (kj), for several sets of initial parameters according to a simple 
batch processing model (see Appendix B). Each panel shows the effect of a change in one of the driving variables 
(lower curve) relative to the amount of energy accumulated for a default parameter set (same upper curve in 
each panel) for which the hypothetical lizard was active for 12 hr at T b - 36*C, activity scope - 1.5, and was 
inactive for 12 hr at T fc - 32*C (see Appendix B for other defaults). PANEL A: the lower curve is for a lower 
inactivity T fc by 6*C to 26*C PANEL B: the lower curve is for a lower activity f » by 4-C to 32*C PANEL C: 
the lower curve shows a 10% increase in activity scope to 1.65. PANEL D: the lower curve shows a~ 10% increase 
in the cost of digestion to 0.85 (see "constant" in eq. B 1 1 Appendix B). PANEL E: the lower curve shows a 10% 
decrease in food availability to 90% of the ad lib, feeding rate. PANEL F: the lower curve shows a 10% decrease 
in food energy density to 23,040 J/g dry wt 



ably an underestimate of that exhibited by 
some lizards (Waldschmidt et ai, 1987). 
According to the model, higher scopes 
quickly lead to zero or negative net accu- 
mulation (not shown), which is consistent 
with field observations of reproductively 
active male lizards (Anderson and Karasov, 
1988; Grant and Dunham, 1990). 



Panel (D) shows that an across the board 
increase in the cost of digestion by 1 0% (e.g., 
if consumed food were more difficult to 
digest, or if our guess for the "constant" in 
eq. Bl, Appendix B were too low) effects a 
10% decrease in the rate of energy accu- 
mulation. Panels (E) and (F) show that 
decreases in estimated food availability or 
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food energy density, respectively, also dra- 
matically affect the energy budget. 

The extreme sensitivity of our energy 
budget model to variation in input param- 
eters as shown in Figure 9 angues strongly 
for directly estimating these parameters 
under field conditions for the particular liz- 
ards of interest. But more importantly, these 
model parameters describe important phys- 
iological and behavioral characters of the 
study lizards and may provide useful bases 
for comparison among populations in dif- 
ferent thermal and food resource environ- 
ments. Formulating appropriately resolved 
batch processing models would more real- 
istically mimic the phenotypic processes of 
energy accumulation by individuals in dif- 
ferent populations. This would provide nec- 
essary optimality criteria to generate 
hypotheses for the evolution of thermal 
physiology and thermoregulatory behavior 
subject to environmental thermal and 
resource availability constraints. We 
develop this idea further in the next section. 

Behavioral optimization 

A substantial short-coming of our present 
macroclimatic constraint/energy budget 
model is the simplistic set of responses 
allowed by the hypothetical ectotherm to 
variation in simulated resource and thermal 
environmental constraints (see Fig. 4); But, 
according to the schematic in Figure l) 
organismal fitness should involve maxi- 
mizing the expected sum of current and 
future fecundity. Toward this end, our model 
must be modified to include explicit pre- 
dictions about how lizards should optimally 
acquire and process resources on a daily 
basis (short term foraging strategy) and opti- 
mally allocate assimilated resource over a 
breeding season and/or lifetime (longer term 
life history strategy). The problem is one of 
reducing a diverse set of environmental 
constraints and risks associated with pos- 
sible "actions" into a single currency, e.g., 
fitness, and formulating an amalgamated 
cost/benefit function that determines the 
time and energy allocation phenotype by the 
hypothetically optimizing organism 
(McQccry, 1978; McNamara and Houston, 
1986; Dunham et ai, 1989). Numerous 
studies have used diverse optimization 



techniques to generate hypotheses about life * 
historical time and energy allocation includ- 
ing exact solutions using optimal control 
(King and Roughgarden, 1982; Iwasa and 
Roughgarden, 1984; Werner, 1986; Horn, 
1988), numerical approximations such as 
the quasi-Newton method (Taylor and 
Gabriel, 1992), or stochastic dynamic pro- 
gramming (McNamara and Houston, 1986* 
Mangel and Clark, 1 988; Horn et al % 1 990)! 

In this section we describe a simple model 
of optimal time allocation by an hypothet- 
ical lizard maximizing its long term survival 
and energy budget size in order to demon- 
strate the utility and ease of inclusion of 
optimality models in our macroclimatic 
constraint/energy budget modeling. We use 
the methods of stochastic dynamic pro- 
gramming since according to Horn et ai 
(1990) these methods are heavily founded 
on ecological specifics of the study system 
and return explicit and testable predictions. 
In our model, our hypothetical lizard may 
be either active or inactive depending on its 
body energy reserves, resource and thermal 
environmental characters, and risk of mor- 
tality due to foraging. As organized by Man- 
gel and Clark (1988), our model has. five 
essential components. First, the "s\ate 
space" is the energy state, x, of the lizard at 
time t. Second, the "constraint set" we use 
is an interval constraint, i.e., x must be in 
[0, max x]. Third, the "strategy set" relates 
simply that the lizard must be either active 
and foraging or inactive (at predetermined 
activity and inactivity T b s, respectively). 
Fourth, the state dynamics relate to how the 
energy state and probability of mortality are 
affected by being active vs. inactive. In our 
model, the lizard's energy state is decre- 
mented each time interval due to its resting 
metabolism (as a function of its T b ), and if 
it is active and foraging, the costs of these 
are also deducted. The lizard's energy state 
can only be increased by being active and 
foraging; however, predators may also be 
encountered during activity, but not inac- 
tivity. Food and predator encounters are 
expressed as probabilities which motivates 
the term "stochastic" to describe this type 
of state variable model. Lastly, the "opti- 
mization criterion" in our model is to choose 
to be active or inactive for each x and t 
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depending on which choice simultaneously 
maximizes the probability of surviving as 
long as possible while attaining the greatest 
energy state after a long time. We present 
the equations and pseudocode for this sim- 
ulation in Appendix C 

Figure 10 illustrates two representative 
response surfaces from our model of opti- 
mal time , allocation (see Appendix C for 
parameter inputs for these simulations). 
Energy state and time are plotted in the hor- 



izontal plane, and since there are only two 
possible choices in the "strategy set*' (active 
or inactive), we plot the predicted fitness 
difference (F. rtive (x, t, T) - F iMCtive (x, t, T)) 
along the Z-axis for each x and t. Positive 
Z-values indicate that activity is optimal, 
and negative Z-values (shaded portions in 
Fig. 10) indicate that inactivity is optimal. 
As can be seen in both response surfaces, 
the model predicts that lizards with low to 
moderate energy states should be active, 
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whereas energetically endowed lizards 
should be inactive to maximize their long 
term fitness. 

The difference between response surfaces 
in Figure 10 is due to an increase in food 
availability (the parameter 0 in Appendix 
C) for which the fitness difference for being 
active increases greatly (note different scales 
of Z), but only for low energy states. We 
have conducted detailed analyses in other 
regions of this six-dimensional input 
parameter space, and the qualitative differ- 
ence between response surfaces in Figure 10 
is typical of variation in other input param- 
eters. In general, there are often large dif- 
ferences in the predicted fitness values as 
well as some variation in where the cutoff 
for activity occurs among states in x; how- 
ever, all simulations were consistent with 
the result that energy state and activity level 
should vary inversely. Additionally, we have 
simulated predation during inactivity (such 
as by fossorially foraging snakes), and we 
have examined the effect of a risk of pre- 
dation during any change from activity to 
inactivity (such as by wait-ambush snake 
predators near crevice openings [we thank 
A. Cannon for this model]), all of which 
support the qualitative prediction in Figure 
10 that activity rates should be lowest for 
well fed lizards. 

According to results from a field resource 
supplementation experinfent, in which 
marked free-ranging Sceloporus merriami 
lizards were individually fed individual 
crickets, the activity rates of fed females were 
significantly lower than werie the rates for 
control unfed females (Fig. 10, Grant and 
Contreras-Arquieta, unpublished data). 
Thus, female S. merriami appear to con- 
form to the predictions of our optimal time 
allocation model. Note however, that the 
activity rates for fed and unfed males did 
not differ (there were no differences in the 
amount of food supplemented to treatment 
females and males). We interpret this dif- 
ference to be indicative of the difference in 
reproductive strategies between the sexes for 
these highly social lizards. For male S. mer- 
riami, reproduction depends greatly on ter- 
ritorial acquisition and exclusion of com- 
peting males, as well as on spatial overlap, 



encounters, and courtships with neighbor- 
ing females (Ruby and Dunham, 1 987). In 
contrast for females, the amount of her 
energy reserves will greatly determine her 
reproductive Output (Dunham, 1981) which 
is consistent- with an underlying assumption 
of the model that energy and predation con- 
straints are sufficient to drive behavior. We 
suspect that the lack of concordance between 
the model predictions and male responses 
is because the assumptions of the model 
poorly described the fitness criteria for male 
reproduction which may be greatly affected 
by operative social/demographic con- 
straints on male reproductive activity (Dun- 
ham al, 1989). 

These results underscore the necessity of 
obtaining a detailed understanding of the 
natural history of the study system prior or 
parallel to any modeling attempt. But more 
importantly, constraints from the social/ 
demographic environment dictated by 
requirements of mate acquisition in polyg- 
ynous social systems can profoundly affect 
adult male behavior and must be included 
in behavioral optimization models. More 
advanced models must recognize that dif- 
ferent sexes and ages may have different 
optimization criteria, and simultaneously 
include all relevant factors impinging on the 
daily behavioral and seasonal resource allo- 
cation phenotype from the biophysical, 
resource, social/demographic and predation 
environments (Dunham et ai t 1989). 

Concluding Remarks * 

The macroclimate constraint/energy 
budget model we describe relies on basic 
physical characters of microclimates and 
well-known biophysical and bioenergetic 
characters of lizards to delineate environ- 
mental constraints on individual activity, 
food consumption, and food processing 
capacity. We stress that this modeling can 
only predict the size of the energy budget 
available for allocation to growth, storage 
and reproduction in a life history and not 
how that resource should be allocated to 
maximize fitness (see Figs. 1 and 3). There- 
fore, models of this type should be used as 
a "front-end" to provide the bounding envi- 
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ronmental constraints on individual ecto- 
therm energy budgets within which future 
models of optimal resource acquisition and 
allocation must operate. We advocate as 
have others (e.g., DeAngelis et aL, 1979; 
Huston et aL, 1988) that such "individual- 
based* 1 models are an important tool to 
explore potential population and commu- 
nity level responses to macroclimatic vari- 
ation. These formulations have the advan- 
tage of being extremely system-specific and 
thereby generate testable predictions about 
the mechanistic effects of interacting envi- 
ronmental constraints on, and the adaptive 
significance of, variation in individual life 
history characters (Dunham et aL, 1989). 

According to many scenarios of climate 
change (COHMAP, 1988) the magnitudes 
of predicted changes in the driving macro- 
climatic characters (e.g., temperature, 
humidity, and cloud cover) are similar to 
presently existing differences among sites 
separated by only a few hundred meters in 
elevation or a few degrees in latitude. There- 
fore, results of interdisciplinary studies of 
present day patterns of covariation in envi- 
ronment and life history among geograph- 
ically distinct populations may be extrap- 
olated to generate hypotheses about 
individual and population level responses 
to climate change in a given location. How- 
ever, such an extrapolation cannot be made 
without a thorough knowledge of the study 
organism and of the complexities of the 
interactions between individuals and their 
total environment (biophysical, resource, 
social/demographic and predation, Dun- 
ham et aL, 1 989). An important lesson from 
numerous among-site comparisons of life 
history and environmental covariation (e.g., 
Sceloporus merriami. Grant and Dunham, 
1990) is that individual responses to inter- 
acting environmental gradients may involve 
significant age, sex, site, and season inter- 
actions among life history characters. 

The correlation between observed and 
predicted energy budget sizes among pop- 
ulations of Sceloporus undulatus (Fig. 6) 
suggests that our modeling efforts are on 
track; however, the limited resolution and 
relatively short duration of most of the com- 
parable field studies as well as the lack of 



data on geographical and seasonal variation 
in physiological ecology are significant 
impediments to modeling advances. Prog- 
ress in this area requires a tremendous 
increase in the resolution, duration and 
number of population ecology studies in 
order to use broad intra- or interspecific 
comparisons to generate and assess hypoth- 
eses concerning the evolution of life histo- 
ries as well as population responses to cli- 
mate change. These studies must be 
interdisciplinary and include detailed phys- 
iological and biophysical modeling as well 
as estimates of constraints from the resource, 
predation and social/demographic environ- 
ments. Such studies must also be of suffi- 
cient duration to elucidate the temporal 
dynamics of individual and environmental 
covariation (Nichols et aL, 1976) since crit- 
ical life history characters such as growth, 
survival and reproduction tend not to be 
constant (e.g., Dunham, 1978, 1981) and in 
fact may exhibit adaptive plasticity (Caswell, 
1983; Stearns, 1986; Newman, 1988). 
Finally, studies must include assessments of 
the genetic bases of individual responses to 
environment in order to model long-term 
evolutionary effects of climate change. An 
important aspect of using interpretations of 
present day geographic covariation in life 
history and environment (for which envi- 
ronments have undoubtedly differed for a 
very long time) to generate predictions about 
the . effects of climate change within a loca- 
tion is to determine the time frame over 
which macroclimatic characters vary rela- 
tive to organismal evolutionary response 
times. Controlled reciprocal transplant and 
common garden experiments are ideally 
Suited to detect genetic effects and should 
be seen as a research priority despite their 
extremely time- and labor-intensive meth- 
ods. 

Clearly, long-term interdisciplinary data- 
bases will be crucial in the development and 
testing of advanced models that link mac- 
roclimate and individual autecology to pop- 
ulation and community ecology, and we 
suspect that the data needed for testing these 
models are likely to exceed the capabilities 
of any individual laboratory in terms of 
labor, logistics, and required degree of spe- 



174 



B. W. Grant and W. P. Porter 



cialization. An organized collaboration 
among research specialists is likely to be far 
more productive than any single generalist 
given the same resources. Therefore, we urge 
in this era of evaporating funding that inter- 
disciplinary research should increase in pri- 
ority as an important means to resist the 
potentially divisive effects of increased 
competition for ever-scarcer research sup^ 
port. 
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Appendix A. 

Equations used in the ectotherm 
physiological simulation . 

Hourly Metabolic Expenditure, hQ mct (J/hr): The 
metabolic rate of an active lizard is assumed to be a 
constant fraction times its resting rate ("activity scope, 
Bennett and Dawson, 1976; Congdon et aL. 1982). 
Resting metabolic rates for this simulation are obtained 
from a regression in Bennett and Dawson (1976): 

Qmet (J/hr) - 20.2- 10^°" - 1 "'»■ . . . 

•lizard mass °"*activity scope (eq. Al) 

The coefficient preceding the first exponent assumes' 
that the lizards are insectivorous (see Bennett and Daw- 
son, 1976), and a value of 1.5 was selected for the 
"activity scope" which is appropriate for active igua- 
nids (Congdon et aL, 1982). 

Hourly Resource Accumulation, »Qavail (J/hr): 
Maximal food consumption and assimilation rates were 
determined for a lizard as a function of T„ by estimating 
the steady state rates of ingestion and excretion (J/day) 
by lizards held at constant cage temperature and in 
either an ad lib, feeding regime or a regime at 50% of 
ad lib. (data and regressions from Sceloporus undulatus, 
see Waldschmidt, 1984; Waldschmidt et aL 1986, 
1987). These rates were then converted to hourly rates 
by dividing by 24 nr. Note that if the lizard is not 
active during a particular hour, then the value for Qavail 
is zero for that hour. The regression for the digestive 
efficiency as a AT*) was obtained from arcsine square- 
root transformed values since these were proportions 
(Sokal and Rohlf, 1981). The generating equations are: 

max consumption rate (mg wet wt/gday) 
- -328.7 + 22.82-T* - 0.320 T» l (eq. A2) 
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actual consumption rate (mg wet wt/hr) 
■ max consumption rate . . . 
fraction of ad Ub. feeding- , . . 
•lizard mass/24 hr 

arcsin digestive efficiency 
-85.34 - 0.50 T b + ... 



(eq. A3) 



f 0 

[o, 



if frac. of ad lib. » i.O 
000074 . T b » if frac. of ad lib. < 1.0 



(eq. A4) 

digestive efficiency 

- [sin[(T/l80) arcsine digestive efficiency]] 3 (eq. A5) 
Qavail (i/hr) « actual consumption rate- 

(1 - food percent water)- . . . 
•digestive efficiency . . . 
* food energy density (cq. A6) 

In eq. A3 the "fraction ofW lib. feeding" ranges 
. between 1 .0 for ad Ub. feeding to 0 for no feeding In 
cq. A6. we assigned values of 0.75 for "food percent 
water* and 25,600 J/g dry weight for "food energy 
density" according to Waldschmidt et at. ( 1 986). 

Appendix B. 

Equations used in the batch 
processing simulation 

Our batch processing resource assimilation model 
uses the specific digestive performance relations derived 
for Sceloporus undulatus (see Appendix A) to find the 
amount of food ingested and assimilated per hour for 
a user defined 24 hour cycle ofT b and activity, and 
then simply sums the total assimilated (which was 
defined earlier as the discretionary resource, »DER) 
over all of the hours during a 10 day period. 

In addition to assumptions mentioned in the text 
(/.*., step functions in activity rate and T b selection in 
time, and food consumption constrained only by gut 
volume), we must also assume some relation for the 
cost of digestion, or the specific dynamic action (= 
Qsda). According to Waldschmidt et al (1987) (L, 
in S. undulatus varies with the amount of food ihthe 
gut as well as with T b , and can be as great as the resting 
metabolic rate </.*.. 1.0 Q MCT ) for lizards with full guts. 
In the absence of any quantitative basis, we assumed 
that Qsda co ; varied linearly with Q MET (eq. Al), and 
was proportional to the ratio of the amount of food 
presently in the gut divided by the maximal gut capac- 
ity (found by differentiating eq. A2 to get the optimal 
T b or 35.7 s , and using this value in eq. A2). Thus we 
assumed. 



formed a polynomial regression of lab data on food 
passage rate for 5. undulatus in Waldschmidt M984 
R* - 0.987): v ' 

passage rale (days) = 102.891 - 8.68I2T„ + 
+ 0.251568-TV - ... 
- 0.00247528-T b > (eq. B2) 
and we assumed that the hourly passage rate was sim- 
ply '/ J4 of the daily rate. Thus, the hourly assimilation 
rate is given by: 

assimilation rate 

„ digestive efficiency (from eq. A5) 
passage rate* 24 hours 



(cq. B3) 



Below we present the pseudocode for this model and 
we encourage readers to construct their own models 
using other versions of the generating equations specific 
to other taxa of interest: 

begin simulation 
begin day 
for each hour of day do 
-if active then get activity tb and metabolic 
scope and fill gut volume, V, to capacity 
(eq. Al); 

else if inactive gel inactivity tb and set met- 
abolic scope equal to 1.0; 

—calculate the hourly assimilation rate (eq. B3); 

—calculate the volume of food digested during 
the hour, d V, from the product of the assim- 
ilation rate and the volume of food in the 
gut, V; 

—calculate the amount of energy in d V, m Q AVAIL 
(J/hr), using Q AVA1L « dV-(l - food percent 
water) food energy density (modified from 
eq. A6, and see Appendix A for coefficients) 

—calculate the metabolic rate, Qmct (eq. A I), 
and the cost of digestion, Q SOA (eq. Bl ); 

—increment the energy budget, DER, by + 

Qavail *~ Qmct - Qsda* 

—decrement gut volume by dV; 
end of hour and increment to next hour ; 
end of day and increment to next day; 
end of simulation; 
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Qsda ~ constant (food volume in gut + . . . 

^ max capacity) * Q MET (eq. Bl) 

We chose a value for the constant in this equation 

/ , /u^V 5 Wilhin lhc ran « c 0 5 - 1 * 0 for A undu- 
latus (Waldschmidt et al.. 1987, see their Table V, p. 
573). We emphasize that data are critically lacking on 
the cost of digestion in lizards. 
To estimate the hourly assimilation rate we per- 



Equations used in the stochastic dynamic 
model of optimal time allocation to 
activity or inactivity by a lizard 

This section describes the specifics of our particular 
model and omits most of the generalities of the meth- 
ods. For readers who are interested in a general descrip- 
tion of the methods of stochastic dynamic program- 
ming, we highly recommend an excellent book by 
Mangel and Clark (1988). In our model, the lizard's 
energy state is decremented each time interval due to 
its total metabolic cost of being active, Q MET {aciive), 
or inactive, Q MCT <inactive), according to equation A I 
for specified activity and inactivity T b s. Further, the 
lizard's energy state can only be increased by being 
active, during which food of value a is encountered 
with probability 0, and predators are encountered with 
probability r (which instantly reduces state and fitness 
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to zero). Thus, if the lizard is active, one of three out- 
comes is possible for its state at time t + I (eq. CI): 

'x(l) ~ Q«ET<active) + a 

w/ probability 0 (1 - V) 
x(t + I) - x(t) - Q MCT <inactive) 

w/ probability (1 - - r) 
0 w/ probability r 

but, if it is inactive, only one outcome is possible: 

x(t + 1) - x(t) - Q MCT <inactive) (eq. C2) 

To encapsulate the optimality criterion, we form the 
dynamic programming equation to find F^^(x t t, T) 
and F^^x, t, T), which are the fitness functions for 
active or inactive lizards, respectively, with state x over 
the interval t to T, and assuming behavior is optimal 
for all times thereafter, i.e., between t + I and T, (see 
Mangel and Clark, 1988, for a highly readable presen- 
tation of the general formulation of the dynamic pro- 
gramming equation) (eq. C3): 

F«iJx, t, T) - [(I - O [0F ( (x[t + I]) + (1 - 0) 
•..••F a (x[t+ 
Fu^iJx, i, T) - F,(x[l + 1]) 

Where, F,(x(t + i]) and F,(x[t + 1]) « the fitness of 
the lizard with state x(t + 1) in the next time interval 
assuming that it was active and successfully found food, 
or was active and did not find food, respectively, 
according to eq. CI. And, F,(xtl + 1J) « the fitness of 
the lizard with state x(t + 1) in the next lime interval 
assuming that it was inactive, according to eq. C2, 
Whether the lizard should be active or not depends on 
the larger of the two fitnesses, F^x, t, T) vs. F^^x, 
i.T). 

Note that these are recursive equations where the 
fitnesses for choices at time t depend on the fitnesses 
at time t + I, thus, these equations can be solved by 
specifying the ending state vs. fitness distribution (also 



called the "terminal fitness function") and iterating 
backwards in time. For our simulation we chose a lin- 
ear terminal fitness function (F[x. T, TJ - x/maximum 
x) as one way to meet our "optimization criterion" that 
the lizard should seek the greatest long-term energy 
state possible during its foraging. 

For the simulation in Figure 10, we choose a 10 min 
interval as our unit of time and 1 J as our unit of energy 
state x, ranging from 0 to 60 J. These selections were 
constrained by our computing capacity; however, we 
note that much levity exists in the appropriate selection 
of time and state resolution since we are interested in 
the qualitative pattern of response (i.e., should the liz- 
ard be active or inactive) and not the absolute fitness 
predictions. We selected activity and inactivity T*s of 
26* and 36*C an activity scope of 2.0 (1 .5 for activity 
+ 0.5 for food processing), a value of 0.05 for predation 
risk r, and a value of 20 J for net food energy content 
a. Lastly, to show the effect on the predicted response 
surface of increased food availability, we used food 
encounter probabilities, 0s, of 0.75 (Fig. 10, upper left) 
and 0.50 (Fig. 10, upper right). 

Below is the pseudocode for this simulation: 

begin simulation 
initialize linear terminal fitness array, 
set fitness for state x [from 0 to max x] 
- x/max x 

initialize the time counter to its maximal value, T 
repeat iterating backwards in lime 
for each energy state from 0 to max x do 
compute the changes in energy state x due to all 
possible outcomes for activity and inactiv- 
ity (eq. CI andeq. C2) 
compute the fitness functions to find the expected 
fitness for choosing to be active or inactive 
(cq.C3) 

if F^ „(x, U T) > Fi^c^x, t, T) then be active, 
otherwise be inactive for x and t 
decrement the time counter 
return to repeat again and quit after a long time into 
the past 

end simulation. ; 
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